Superemitters in Hybrid Photonic Systems: A Simple Lumping Rule for 
the Local Density of Optical States and its Break-Down at the Unitary Limit 
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We theoretically investigate how the enhancement of the radiative decay rate of a spontaneous emitter pro- 
vided by coupling to an optical antenna is modified when this "superemitter" is introduced into a complex 
photonic environment that provides an enhanced local density of optical states (LDOS) itself, such as a micro- 
cavity. We show that photonic environments with increased LDOS further boost the performance of antennas 
that scatter weakly, i.e. that are far from the unitary limit, for which a simple multiplicative LDOS lumping rule 
holds. In contrast, enhancements provided by antennas close to the unitary limit, i.e. antennas close to the limit 
of maximally possible scattering strength, are strongly reduced by an enhanced LDOS of the environment. Thus, 
we identify multiple scattering in hybrid photonic systems as a powerful mechanism for LDOS engineering. 



Optics encompasses the most fascinating part of the elec- 
tromagnetic spectrum, due to its energetic overlap with elec- 
tronic transitions in matter. Nano-photonics aims at con- 
trolling such transitions by molding light at sub-wavelength 
scales. Purcell first predicted that resonators modify the ra- 
diative lifetime of spontaneous emitters, a property held to be 
intrinsic to the source until then [1]. Modern literature dis- 
cusses the Purcell effect in terms of the local density of op- 
tical states (LDOS), a fundamental quantity governing spon- 
taneous emission, thermal radiation, and vacuum forces . 
Two tools have emerged to shape the LDOS: on the one hand, 
interfaces Jit], photonic crystals fl-El and dielectric micro- 
cavities 10] modulate the LDOS on length scales of order A/2 
via interference. On the other hand, optical antennas 10 ^ 
employ intrinsic plasmonic material resonances to enhance 
the LDOS on A/20 length scales. Optical antennas are so 
small compared to the wavelength A that a source-antenna 
ensemble essentially radiates as a dipole. This similarity to 
a bare molecule, but with much higher radiative rate, has 
prompted the term "superemitter" 111 Oil to refer to this source- 
antenna entity. The availability of photonic building blocks on 
such different length scales raises the exciting idea of integrat- 
ing deep- sub wavelength superemitters in much larger dielec- 
tric structures to obtain a combined advantage , e.g . by placing 



a nano-antenna inside a micro-resonator fllli 11211 . In view of 



these recent developments it is imperative to understand how 
the LDOS of such hybrid systems emerges from that of the 
separate entities. Since nano-optic devices can be interpreted 
as lumped optical elements HJ, , where the LDOS acts as 
an impedance for spontaneous emitters, much like the resis- 
tivity experienced by a classical AC current source due to the 
fact that it radiates energy into the far field, some simple cir- 
cuit rule might be hoped for, which lumps the LDOS provided 
individually by the photonic building blocks. 

This Letter theoretically investigates how the decay rate en- 
hancement provided by a superemitter varies when it is posi- 
tioned within a larger photonic system, i.e., how the LDOS 
lumps. We consider two archetypical examples: First, we in- 
vestigate a superemitter coupled to high-Q resonances. Sec- 
ond, we consider a superemitter in front of a perfect mirror. 
Our analytic, yet fully electrodynamic model yields a simple 



multiplicative LDOS lumping rule for moderate antenna fac- 
tors. For antennas operating close to the unitary limit, how- 
ever, this simple lumping rule breaks down and the total en- 
hancement becomes proportional to the inverse LDOS of the 
background system. This insight paves the way for engineer- 
ing the LDOS by exploiting multiple scattering in hybrid pho- 
tonic systems. 

The radiative decay rate enhancement of a spontaneous 
quantum emitter can be calculated via the power required 
to sustain a constant classical point current j = p Q that 
loses energy by radiation |0]. This power equals the cycle 
averaged work per unit time done by the source's electric 
field on its own dipole moment p . The electric field gen- 
erated at position r due to a source p at vq is calculated 
via the electric Green function G(r,ro) of the respective 
system. This yields the power required to drive the source 
P = ^oupQlmG(ro,ro)p . For lossy environments, this 
expression describes the total decay rate, i.e., radiation plus 
quenching induced by the environment. We use the term 
LDOS to refer to ImG projected on a unit vector along p Q , 
i.e. to the decay rate modification of a molecular dipole p 
at position r o . Every strategy to boost decay rates via a pho- 
tonic structure, be it large and dielectric, or a nano-antenna, 
represents a modification of ImG. The small size and dipolar 
nature of a nano-antenna however suggest to interpret its rate 
enhancement rather as an enhancement of the dipole moment 
p instead of a change in ImG fioh . The simplest optical an- 
tenna is just a particle with polarizability tensor ol(uj) (8|]. At 
distance d in the near field of an emitter, the particle acquires 
a large dipole moment ex oc/d 3 . The total dipole moment of 
the emitter-particle assembly p = [1 + oc/d 3 )po can exceed 
by far the source's dipole moment p Q , rationalizing the term 
'superemitter'. If d <C A, the power radiated by the superemit- 
ter itH 



\u |po| 2 LDOS B (r ) x A 



(1) 



exceeds that of the bare source by the antenna factor A — 



|l + a/d 3 



\a\ /d 6 . The rate enhancement provided by 



the embedding background system at the location of the su- 
peremitter is described by LDOSb- Mie calculations show that 
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A accurately describes antenna particles up to 60 nm in diam- 
eter in vacuum lfl6ll . Equation (Q} suggests that a superemitter 
can be used as replacement of a bare molecule to probe the 
LDOS of a larger photonic system, since in this reasoning a 
simple product rule lumps the enhancements provided by the 
antenna A and by the photonic environment's LDOSr. 

We analyze the hypothesis of a simple lumping rule in an 
analytic electrodynamic point scattering theory, which is ex- 
act to all multiple scattering orders, with the sole assump- 
tion that the scatterers that constitute the antenna are well ap- 
proximated as point dipoles. The scatterers 1 ... N acquire 
dipole moments p ± . . .p N in proportion to their polarizabil- 
ities cti . . . ocjsf, and the electric fields E(r\) . . . E(tn) at 
their locations r„,, according to the linear self-consistent set 
of equations to, \lm 



Superemitter ! Hybrid: Superemitter + Background 



Pn <*r. 



E x 



(2) 



By construction, the antenna described by the cx n is explicitly 
separated from the background system that it probes, which is 
quantified by its Green function Gb . For a consistent theory, 
three facts need to be taken into account. First, the particle 
polarizability directly depends on the background via fll4l1 



a n0 ~ G B (r n ,r n ), 



(3) 



where ol Uj o is the electrostatic polarizability. Note that 
isotropic particles can acquire an anisotropy due to the 
anisotropy in radiation damping given by ImGs and reso- 
nance energy shift due to ReG# in complex photonic sys- 
tems Ill9ll . Second, the source in our model is a single dipole 
p at r (with a = 0) so E- m (r n ) = G B (r n , r )p . Third, 
the total decay rate of the source in the full system is found 
via the cycle- averaged work done by the total electric field on 
the source p yfl . Therefore, using a source of unit strength 
p , the LDOS equals 

pi • ImG B (r , ro ) • P + Im Yl p ® ' G #( r o> r n)-P n - ( 4 ) 

n>l 

The first term is the LDOSb provided by the background sys- 
tem itself in absence of the antenna, while the second term 
is the contribution from the antenna. We calculate the LDOS 
of the hybrid system by evaluating Eq. ©, after solving the 
37V linear equations in Eq. © for the induced dipole moments 
p x . . . p N . We use the exact Green function for a sphere lEoll 
and a planar interface Q |2l]] to evaluate how a superemit- 
ter probes the two canonical cases of micro-cavity resonances 
and non-resonant interfaces yfl, respectively. 

As a superemitter we consider a source in the gap of a 
dimer antenna with its dipole p along the symmetry axis [see 
FigOIa)]. The Ag spheres forming the dimer, each of diam- 
eter 40 nm, have a center-center distance of 60 nm [12211 . The 
dashed line in Fig.QJb) shows Ima of a single antenna particle 
in vacuum. The solid line is the decay rate enhancement of the 
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FIG. 1. (a) Left: Sketch of superemitter formed by two silver spheres. 
The fluorescent source is located between the two particles with its 
dipole moment along the symmetry axis. Right: Hybrid photonic 
system comprised of superemitter embedded in background system 
formed by dielectric sphere. Red and blue colors illustrate fields of a 
whispering gallery mode, (b) Dashed line: polarizability of a single 
antenna particle in vacuum. Solid line: antenna enhancement factor 
for superemitter sketched in (a) in vacuum, (c) Dashed line: Extinc- 
tion efficiency of dielectric sphere showing narrow Mie-resonances. 
Solid line: Purcell factor of the Mie-sphere 50 nm from its surface. 



source in the antenna gap with the antenna in vacuum. The an- 
tenna factor reaches about 1200, in good agreement with full 
multipole calculations 1I23I1 . which justifies the use of a dipole 
model. The dimer resonance is red- shifted and broadened by 
super-radiant damping compared to the single particle due to 
longitudinal symmetric dipolar plasmon hybridization ll24ll . 

As a background system we now consider a glass sphere 
(n=1.5) with diameter 2.4 um, supporting whispering gallery 
modes, illustrated in Fig. QIa). The dashed line in Fig. Q2c) 
is the sphere's extinction efficiency J25I1 . Characteristic reso- 
nances are clearly visible as sharp peaks with Q w 10 2 . The 
Purcell factor for a radially oriented source 50 nm from the 
sphere surface [solid line in Fig. UXc)] reaches moderate val- 
ues around 10. Only every second peak in extinction appears 
in the Purcell enhancement, which reflects the field orientation 
of the modes according to the common TE/TM type classifica- 
tion Il26ll . We now place the superemitter with the center of the 
closest antenna particle 50 nm from the sphere's surface and 
the symmetry axis pointing radially outwards [see Fig. QIa)]. 
The solid line in Fig. [2a) is the decay rate enhancement expe- 
rienced by the emitter embedded in the dimer antenna which 
in turn is located next to the Mie-sphere, i.e. the rate in the 
hybrid system normalized to the rate of the bare source in vac- 
uum. While the overall shape of the enhancement provided by 
the antenna in vacuum [Fig.|2a), dashed] is still visible, sharp 
features appear at five spectral positions coinciding with the 
sphere's whispering gallery modes [cf. Fig. Q2c)]. To illus- 
trate the effect of the background system on the superemit- 
ter we now normalize the decay rate enhancement provided 
by the hybrid system [solid line in Fig. EJa)] to the enhance- 
ment provided by the bare antenna in vacuum [dashed line in 
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FIG. 2. (a) Solid line (hybrid system): Decay rate of source in the 
gap of the superemitter next to the Mie- sphere normalized to rate of 
source in vacuum. Dashed line (superemitter in vacuum): Decay rate 
of source in the gap of the superemitter in vacuum normalized to rate 
of source in vacuum. At the Mie-resonances, the hybrid LDOS en- 
hancement is drastically modified compared to that of the antenna 
in vacuum, (b) Solid line (hybrid system): Decay rate enhancement 
in hybrid system [solid line in (a)], normalized to the rate enhance- 
ment in the superemitter in vacuum [dashed line in (a)]. Off antenna 
resonance, the superemitter benefits from the LDOS enhancement 
offered by the sphere (dashed line), while on antenna resonance the 
enhancement is suppressed by the inverse of the sphere's LDOS (dot- 
ted line). 



Fig.Oa)] and plot it as the solid line in Fig. Ob). The sharp 
enhancements located in the wings of the antenna resonance 
follow the LDOS of the sphere, denoted by the dashed line in 
Fig. [2b). Therefore, off antenna resonance, at still significant 
antenna factors, we indeed find the anticipated behavior of a 
superemitter, i.e., its already enhanced decay rate compared to 
the bare molecule is further boosted by the presence of a high- 
Q resonance of the background. Furthermore, we note the dis- 
persive features in the enhancement, which swap orientation 
upon crossing the antenna resonance. Surprisingly, however, 
on antenna resonance, the LDOS enhancement is strongly 
suppressed by the Mie-sphere. This LDOS suppression close 
to antenna resonance cannot be explained by a spoiling or de- 
tuning of the cavity by the antenna, since this would only re- 
sult in a shift or absence of a sharp line of extra enhancement 
on top of the bare antenna factor. Also, Waldron's formula 
predicts that our antenna does not significantly shift or spoil 
the micro- sphere resonances due to their large m ode volumes 
compared to a and their low Q-factors |26|,|27|]. Figure Wfc) 
hence implies a spoiling of the antenna by a cavity resonance 
tuned close to the antenna resonance, instead of vice versa. 

In order to understand the truly surprising spoiling of the 
antenna enhancement by a large background LDOS, e.g. of 
a micro-cavity, one needs to interpret Eq. (Q]> correctly, by 
taking the proper radiation damping correction according to 
Eq. © into account. For simplicity we now discuss an an- 
tenna significantly smaller than the wavelength and described 
as a single scatterer with polarizability a [ 2811 . Equation © 
ensures that the optical theorem Ima > LDOSb |a| 2 is 
fulfilled, where equality holds for the case of no material 
loss J29ll . Thereby Eq. © strictly bounds the polarizability 
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FIG. 3. (a) Real (solid) and imaginary (dashed) parts of the ra- 
dial component of the Green function of the Mie-sphere 50 nm from 
sphere's surface. These terms enter the radiation correction to c*Mie. 
(b) Solid line: Radial component of the antenna polarizability aMie 
when located 50 nm from the Mie-sphere. Dash dotted line: a vac 
for same antenna in vacuum. While a vac is limited by the inverse 
vacuum LDOS (dashed line), aMie is bounded by the inverse of the 
sphere's LDOS (dotted line), which leads to the suppression of aMie 
close to the antenna resonance. 



according to 



Im(a) < 



IitiGb 



(5) 



which is a general form of the unitary limit, a fundamental re- 
lation in any scattering theory [f3Qh . In vacuum, this limit is 
well-known for extinction cross section as cr ex t = 47r/clma < 
^:A 2 , a limit reached by an ideal scatterer on resonance, and 
closely approached by any plasmon particle above 20 nm in 
size IIIctI . The unitary limit Eq. © states that a of a strong 
scatterer, and hence the dipole moment it acquires, is propor- 
tional to the inverse LDOS. Intuitively, since a strong scat- 
terer is predominantly damped by radiation, increasing the 
background LDOS further increases its loss and therefore sup- 
presses the scatterers response, which reflects in a reduced po- 
larizability with increased spectral width Bill . 

To quantitatively verify that the unitary limit indeed 
governs the hybrid system's LDOS we evaluate Eq. © 
for a physically small superemitter. The hybrid sys- 
tem's LDOS is then dominated by ReG V ac to read 
ReG va c(^o 5 ri)ImaReGvac(Ti 5 ^o)» which is of order 1/d 6 
(with d = \ri — To |) and precisely yields Eq. (Q} ll32ll . To 
illustrate the effect of the background system on ImaMie of 
a particle located in close vicinity of the sphere we plot in 
Fig- Eta) the radial component of the sphere's Green function 
GMie 50 nm from the surface. Real and imaginary part of Gmiq 
show the typical dispersive and dissipative line-shape of a res- 
onance, respectively. Naturally, ImGMie equals the micro- 
cavity Purcell factor at the source position [cf. Fig. Q2c)]. 
In Fig. [3tb) we plot as the solid line the radial component 
of ImaMie of an antenna located 50 nm from the sphere sur- 
face. Note that the values of Gmiq in Fig.Oa) are the correc- 
tion terms entering Eq. © causing the modification of aMie 
close to the sphere [solid line in Fig. Ob)] with respect to 
a vac in vacuum [dash-dotted line in Fig. [3b)]. Close to an- 
tenna resonance ImaMie is indeed limited by the inverse of 
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FIG. 4. Relative LDOS enhancement for superemitter as a function 
of distance to a near-perfect mirror, (a) Superemitter oriented parallel 
to mirror surface [see inset]. Solid line: Superemitter with source far 
below antenna rensonance [u = 4.0 x 10 15 s~\ cf. Fig.CQb)], traces 
the LDOS of the mirror [dashed line], except at small superemitter- 
mirror separations. A source close to antenna resonance [long dashed 
line, uj = 4.5xl0 15 s _1 ] traces the inverse of the mirror LDOS [dot- 
ted line], except at small separations, (b) Same as (a) for superemitter 
perpendicular to mirror [see inset]. 

the sphere's LDOS [Fig. Ob), dotted line]. Off antenna res- 
onance, the correction to the broad cxmiq due to the narrow 
ReGMie creates characteristic Fano-resonances [33]] . The tran- 
sition from enhancement to inhibition in the lumped LDOS 
between the limits of weak and maximally strong scatter- 
ing is captured by carrying out the correction according to 
Eq. © in the antenna factor A in Eq. (Q}, such that when 
neglecting the real frequency shift the radiated power reads 
P oc | C£q 1 — i LDOSb I LDOSb. For small olq, i.e. in 
the limit of weak scattering, P oc |o:o| 2 LDOSb since weak 
Rayleigh scatterers are constant current sources unaffected by 
the unitary limit ll34ll . Therefore, for weak scatterers a simple 
multiplicative lumping rule for A and LDOSb does hold. In 
the limit of strong scattering P oc LDOSg 1 , since a scatterer 
at the unitary limit is not a constant current source l36ll . 
Accordingly, close to antenna resonance, the enhancement 
curve of the hybrid system [solid line in Fig.[2{b)] follows the 
inverse of the Mie- sphere's LDOS [dotted line in Fig. Ob)]. 

To illustrate that our results are generally valid beyond the 
specific case of high-Q Mie resonances, we examine a near- 
perfect mirror (e=-200) as a background system that modifies 
the LDOS without any resonances [3]. In Fig. Ufa), we plot the 
decay rate of the superemitter in front of the mirror normal- 
ized to the superemitter in vacuum as a function of distance 
to the mirror. The antenna axis is parallel to the mirror [see 
sketch in inset] . Far below antenna resonance, the enhance- 
ment [Fig.|4ja), solid line] follows the mirror LDOS, plotted 
as the dashed line, as expected from the multiplicative LDOS 
lumping rule. Close to antenna resonance [long dashed line], 
the antenna is close to the unitary limit and therefore the to- 
tal rate enhancement follows the inverse of the mirror LDOS, 
indicated by the dotted line. At intermediate frequencies we 
observe a smooth transition between the two cases illustrated 
in Fig. H] (data not shown). For the superemitter oriented per- 
pendicularly to the mirror [Fig. fflb)] we observe the analog 
behavior as for the parallel case. The product lumping rule 
and its cross-over to inverse proportionality are hence generic. 



The framework only breaks down for very small separations 
between superemitter and mirror, where the superemitter size 
is comparable to its distance to its own mirror image. Here 
superemitters also sense gradients in G, which may provide 
a tailorable analogon to recent experiments on mesoscopic 
quantum dots probing multipolar LDOS fl37ll . 

In conclusion, we have examined how the LDOS inside a 
superemitter probes the LDOS of a complex photonic environ- 
ment. Generally for any superemitter with a moderate antenna 
factor the LDOS enhancements of antenna and background 
multiply and a small superemitter can therefore serve as an 
LDOS probe for a large background system, exactly as the 
term suggests ifioll . In surprising contrast, a superemitter with 
an antenna at the unitary limit probes the inverse background 
LDOS, since increasing radiation damping reduces the polar- 
izability of strong scatterers ff30h . Our findings imply that if a 
general lumping rule for optical source impedances fll4ll . i.e. 
LDOS, can be found at all, it must take into account not just 
the bare superemitter LDOS, but also how close the antenna is 
to the unitary limit, which in turn depends on the background 
system. Our insights bear a plethora of exciting prospects 
for en gine ering of LDOS with the toolbox of multiple scat- 
tering 113 80 . One might consider nano-manipulative switching 
of superemitters by moving them with respect to high-Q res- 
onators 113 911 . Furthermore, our work could lead the quest for 
ultra-strong optical antennas towards counterintuitive hybrids 
of nano-antennas embedded in photonic bandgap devices. As 
a further implication, our findings shed new light on attempts 
to use a scatterer as a broadband probe of thermally populated 
modes, since the scatterer acquires spectral features thanks 
to the environment that is to be probed [40]. Notably, hy- 
brid photonic systems composed of superemitters operated off 
resonance together with a micro-resonator might benefit from 
large field enhancements inside the antenna that occur with 
convenient moderate-Q cavities right at the Fano resonances 
in hybrid LDOS, an exciting prospect for single molecule de- 
tection ifiUl . 

This work is part of the research program of the "Stichting 
voor Fundamenteel Onderzoek der Materie (FOM)", which is 
financially supported by the "Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek (NWO)". 
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